A pseudointegrable Andreev billiard 
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A circular Andreev billiard in a uniform magnetic field is studied. It is demonstrated that the 
classical dynamics is pseudointegrable in the same sense as for rational polygonal billiards. The 
relation to a specific polygon, the asymmetric barrier billiard, is discussed. Numerical evidence is 
presented indicating that the Poincare map is typically weak mixing on the invariant sets. This link 
between these different classes of dynamical systems throws some light on the proximity effect in 
chaotic Andreev billiards. 



PACS numbers: 05.45.-a, 74.50.+r 
I. INTRODUCTION 



Billiards have played a prominent role in the under- 
standing of classical and quantum mechanics. In such a 
system, a particle moves freely in a domain with spec- 
ular reflections at the boundary (the angle of reflection 
equals the angle of incidence); see, e.g., JO]. Billiards of 
a different kind are realized in ballistic mesoscopic sam- 
ples connected to a superconductor [@||. The bound- 
ary of such an Andreev billiard Q consists of normal- 
conducting regions with specular reflections and super- 
conducting regions with Andreev reflections, whereby 
electronlike quasiparticles with charge —q, mass to, and 
energy e above the Fermi energy in the normal metal 
are retroreflected as holelike quasiparticles with charge q, 
mass —to, and energy — e. In the absence of a magnetic 
field, retroreflected orbits are self-retracing and therefore 
periodic. The presence of a magnetic field allows for a 
richer spectrum of dynamical behavior |Q] . 

In this paper, we study an interesting Andreev bil- 
liard, the circular Andreev billiard in a uniform magnetic 
field. Our analysis will lead one to the conclusion that 
boundary points separating normal and superconducting 
regions, henceforth called critical points, have the same 
consequences on the classical dynamics as critical cor- 
ners in rational polygonal billiards. In such a polygon, 
all angles otj — mjir/rij between sides are rationally re- 
lated to 7T, where rrij, rij > are relatively prime integers. 
The free motion inside a rational polygon is characterized 
as pseudointegrable Q since it shares some properties of 
integrable systems: (i) the phase space is foliated by two- 
dimensional invariant surfaces |],|7]; (ii) the flow on these 
surfaces is ergodic and not mixing and, in particu- 
lar, not chaotic (see, e.g., ||] for the definition of ergodic 
properties). Yet, in the presence of critical corners with 
nij > 1 the dynamics is more complex: (i) the genus 
of the surfaces is greater than one J5|; (ii) the dynamics 
is not quasiperiodic; and (hi) presumably weak mixing, 
but this is proven only for a special subclass; see, e.g., 
Numerical evidence for weak mixing has been reported 
in 

The paper is organized as follows. In Sec. |l], we intro- 
duce the circular Andreev billiard and derive its Poincare 
map, which describes the collision-to-collision discrete 



dynamics. We show that the dynamics is pseudointe- 
grable, in the same sense as for rational polygons. The 
relation to the asymmetric barrier billiard is discussed in 
Sec. Ill . Section |y| presents numerical evidence that the 
Poincare map is typically weak mixing on the invariant 
sets. In Sec. we draw conclusions and give an outlook. 



II. THE CIRCULAR ANDREEV BILLIARD 



Let us first consider the conventional circular billiard 
in a magnetic field with strength B directed perpendicu- 
lar to the plane. The classical motion of a particle with 
mass to, charge q, and speed v confined inside a circle 
with specular reflections is integrable; see Ref. jl^] and 
references therein. The orbits consist of a series of arcs 
of circles with the Larmor radius R — mv/(qB). With- 
out loss of generality, we scale the radius of the billiard 
and the absolute value of the momentum to unity. If 
R < 1, some orbits form complete circles entirely in- 
side the boundary. Ignoring these complete circles, we 
can specify each orbit by giving the sequence of its po- 
sitions and directions immediately after each impact at 
the boundary. The position on the circular boundary 
is parametrized by the arclength <f> 6 [0,27r). The di- 
rection of the orbit after impact is labeled by the angle 
of reflection a € [0, it], or by the tangential momentum 
p = cos a e [—1,1]. Elementary geometry depicted in 
Fig. |l| gives the Poincare map, i.e., the discrete bounce 
map from the nth to the (n + l)th collision with the 
boundary 



4>n+i = <t>n + u(p n ) (mod 2ir) , 

Pn+1 = Pn , 



with 



u>{p) = 2 arctan 



RVT 



p 



l + Rp 



(1) 
(2) 



(3) 



modulo 2ir restricts the variable to the interval [0,2%). 
The function u{p) is illustrated in Fig. ^. Note the bro- 
ken time-reversal symmetry, cj(— p) =/= —u)(jp). 
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FIG. 1. Part of a typical trajectory (solid arcs of circles) with 
specular (S) and Andreev (A) reflections in the circular Andrecv 
billiard. Dashed lines serve for the construction of the maps (pi)— (H). 




FIG. 2. ui as function of p according to Eq. (]3|) with R > 1 
(solid), R = 1 (dashed), and R < 1 (dotted). The lines u) = n and 
u) = —it are identified. 

The properties of the map (|l|)-(||) are related to the 
continuous-time evolution in a simple manner: conser- 
vation of p corresponds to conservation of angular mo- 
mentum; ergodic motion on an invariant circle p = const 
for irrational u) is related to ergodic motion on a two- 
dimensional invariant torus; families of fixed points for 
rational to correspond to resonant tori foliated by peri- 
odic orbits. 

The situation with a superconducting interface at <j) G 
(0, (3) is illustrated in Fig. |l|. For simplicity, let us first 
assume that the quasiparticles are exact at the Fermi en- 
ergy, i.e. e = 0. An Andreev reflection at the interface is 
then just a change of sign of the tangential momentum 
p; the replacement of an electronlike quasiparticle by a 
holelike quasiparticle with the same energy and vice versa 
can be ignored since the simultaneous change of the sign 
of the charge and the mass does not alter the dynamics. 
Incorporating the change of the sign of p at the interface 
into the Poincare map (0)-(§) gives 

4> n + Lo(pn) (mod 2tt) , (4) 

(5) 



Pn+1 = Pn${4>n) 



with 



The tangential momentum p is no longer a constant of 
motion, but \p\ is. An invariant set \p\ = p > con- 
sists of two circles p = po and p = —po, which are sep- 
arated in phase space (<fi,p). The topology of this situ- 
ation is schematically illustrated in Fig. ||a. An invari- 
ant surface of the continuous-time dynamics is obtained 
from the invariant set by attaching circles as shown in 
the sequences of Figs. |^-||d. The circles represent the 
radial motion which is not contained in the Poincare 
map: after leaving the boundary, the radial coordinate 
decreases until it reaches its minimum; then it increases 
again until it reaches its maximum at the boundary. At 
a normal-conducting boundary point, this path is a full 
loop. Hence, we attach such a circle to each point of the 
set (f> G (/?, 2n) and p G {po, — Po} (solid lines in Fig. ||b). 
At a superconducting boundary point, p changes to — p 
and we have to trace the path once again to obtain a 
full loop. Hence, we connect each point <fi G (0, j3) and 
p = Po with the opposite point p = — pq (dashed lines 
in Fig. ||c) and vice versa. The resulting surface has the 
topology of a two-handled sphere (genus 2) as illustrated 
in Fig. ||d. At the critical points <f> = and tfi = /3 the mo- 
tion is not well defined. Two neighboring orbits hitting 
the boundary of the billiard at different sides of a critical 
point separate from each other by moving along different 
handles of the invariant surface. This is fully analogous 
to the situation near critical corners in rational polygonal 
billiards; cf., e.g., 




FIG. 3. Construction of the invariant surfaces from the invari- 
ant sets. Dotted lines mark the superconducting interface. Thick 
dots mark the critical points. 

Let us now briefly demonstrate that the situation does 
not change qualitatively when the quasiparticles arc not 
exact at the Fermi energy, i.e. e > 0. The electronlikc 
quasiparticle with tangential momentum p e is reflected 
in a holelike quasiparticle with ph ^ — p e - Both parti- 
cles have different Larmor radius, resulting in a different 
to. However, if we redefine the tangential momentum as 
p = ap + b with a = 2c/(p e - Ph), b = -a(p e +ph)/2, 
< c < 1, and also redefine w correspondingly then we 
recover the map (§)-(||). Hence, it is sufficient to consider 
e = as we will do in the following. 



-1 if < (j> < /3 
1 otherwise. 



(0) 
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III. THE ASYMMETRIC BARRIER BILLIARD 

Before analysing the map (|4|)"@ in more detail, we 
discuss an interesting relation to a specific rational polyg- 
onal billiard. A particle with unit mass moves freely in- 
side a polygon consisting of a vertical line of length b 
placed in a rectangle with width L~ + L + and normal- 
ized height 1; see Fig. [|. The symmetric case L~ = L + 
with b = 1/2 is the usual barrier billiar d [|13| - [l5| ; more 
general cases have been considered in |iq-|l8||. Again, the 
trivial energy dependence is scaled away by setting the 
energy to 1/2, or equally the magnitude of the momen- 
tum (pxiPy) to 1. Starting with an initial momentum, 
only a finite number of directions can be achieved during 
time evolution, stemming from the fact that all angles in 
the polygon are rational multiples of n. In phase space, 
the motion takes place on two-dimensional invariant sur- 
faces (Ipzl, |py|) = const. The general formula for the 
genus of such surfaces || gives two due to the critical 
corner at the end of the barrier. 



L 



FIG. 4. Asymmetric barrier billiard, rectangle with a vertical 
line connecting the origin of the coordinate system (x, y) = (0, 0) 
with the point (x,y) = (0, ft). 

It is convenient to consider the barrier billiard as two 
rectangular billiards, one with x > and one with x < 0, 
connected by the passage x — 0, y > b. Suppose for 
a moment that the passage is closed, 6=1. For each 
of the two integrable rectangular billiards we introduce 
action variables I x — \p x \L/tt and I y — \p y \/^ where 
L stands for L~ and L + . The time dependence of the 
angle variables is given by 4> x (t) = u> x t (mod 2ir) and 
4> y {t) = uiyt (mod 2ir) with the frequencies uj x — ir 2 I x /L 2 
and uiy = TT 2 I y . The flow on the torus is ergodic if and 
only if the winding number p — io y juj x — I y L 2 /I x is irra- 
tional. When the passage is open, b < 1, the particle can 
move from one rectangle to the other one. We label the 
rectangles with the sign of x, s = ±1. We then introduce 
the Poincare section x = 0, i.e. we look at the line in con- 
figuration space where s(t) possibly changes. Choosing 
the origin of the angle variables and the barrier length b 
such that the passage is given by < <f) y < (3 = 2n(l — b) 
we get the map 



f>y,n+i = <t>y,n + 2np(s„) (mod 2n) 
s n +i = s„ <!>(</>„) . 



(7) 
(8) 



following sense. Consider an orbit (<fo,Po), (<^i,Pi), • ■ ■ m 
the Poincare map of the Andreev billiard De- 
note the two "frequencies" as lo + = to(\p\) mod 27r and 
uj~ = u>(—\p\) mod 2-7T. The family of orbits parametrized 
by all accessible initial conditions (4>o,Pq) (with fixed ui + , 
u>~) is an invariant set and related to an invariant sur- 
face of the Andreev billiard as described in Sec. [0]. Con- 
sider now an orbit ((f> y ,Q, sq), (4> y ,i, si), • ■ ■ in the map (]?])- 
(|) with 27rp(+l) = uj + and 2?rp(-l) = u~ . This 
can be achieved by constructing a barrier billiard with 
L + = uj + \p x /p y \/(2ir) and L~ — lu~ \p x /p y \/ (2it) where 
(Px-,Py) is the initial momentum. The family of orbits 
parametrized by all accessible initial conditions (4> y ,o, sq) 
(with fixed (p x ,p y ), L + , and L~) corresponds to an in- 
variant surface of the barrier billiard. The two considered 
families of orbits are identical if we identify <p y with </>, 
and s with sign(p). This interesting relation between 
orbits indicates that invariant surfaces in the circular 
Andreev billiard not only have the same topology as in 
the barrier billiard, but also the dynamics on these sur- 
faces (restricted to the chosen Poincare surfaces of sec- 
tion) have typically the same ergodic properties. We will 
discuss this issue in detail in the following section. It is 
to emphasize that for our purpose it is not relevant that 
this kind of equivalence is not complete (we have ignored 
the trivial circles lying entirely inside the billiard bound- 
ary for R < 1) and possibly not bijective (i.e. given a 
family of orbits in the barrier billiard there may be no 
counterpart in the Andreev billiard). A complete corre- 
spondence between barrier and circular Andreev billiard 
is not expected since the former one is symmetric under 
time reversal whereas the latter one is not. 

Note that both maps have four points at which the dy- 
namics is discontinuous, (<fi y ,s) = (0, ±1) and (4> y ,s) = 
([3, ±1). In the Andreev billiard these points are two 
copies of the two critical points; see Fig. ||a. In the barrier 
billiard, the points are four copies of the critical corner. 



IV. DYNAMICS ON THE INVARIANT SETS 



We now discuss the dynamics of the piecewise lin- 
ear, area-preserving map (Q) — (JH|) on the invariant sets. 
Clearly, the dynamics is not chaotic, since all Lyapunov 
exponents are zero. However, weaker ergodic properties, 
such as, mixing, weak mixing, and ergodicity may be 
present. First, ergodicity on the invariant sets follows 
directly from the fact that an orbit of the map (|J)-(|^) 
has a counterpart in the map ((?])- (||) which is typically 
ergodic because the flow on the invariant surfaces in ra- 
tional polygons is ergodic and not mixing [^J. However, 
mixing behavior of the map (Q)-(H) cannot be excluded 
with this reasoning. We do this numerically by iterating 
Eqs. (|J)-(H) and computing the time-averaged autocor- 
relation function (AF) 



Even though the functions 27rp and lo are quite differ- 
ent, the map (0)-(§) is related to the map (|)-(|) in the 



Rz(J) 



%n-\-j 



z n ) 



(4) 



(9) 
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where z stands for p — (p) and <fi — ((f)), respectively. We 
have found that the AF typically does not decay to zero 
which excludes the mixing property. However, the inte- 
grated AF 



1 n 

R,m(") = -Y,\R*U)\ S 

n — ' 



(10) 



3=0 



is found to vanish according to a power law, n~ D2 with 
< £>2 < 1 for large n; D2 is the correlation dimen- 
sion of the spectral measure . The example in Fig. || 
shows that the integrated AF for both p n and <j> n clearly 
obeys the power law with D2 ~ 0.513 and D2 ~ 0.577, 
respectively. We observe qualitatively identical behavior 
for several nontrivial functions f(p n , 4>n) m a more or less 
pronounced way which is consistent with weak mixing as 
maximal ergodic property. 

However, this does not directly imply that the 
continuous-time evolution is also weak mixing on the 
invariant surfaces. For example, the barrier billiard in 
Fig. ^ is not weak mixing since y(t) is a periodic func- 
tion. Nevertheless, since we cannot find such a trivial 
component in the continuous-time evolution of the An- 
dreev billiard, we believe that it is weak mixing, but this 
certainly needs further studies. 




2 4 6 
ln(n) 

FIG. 5. Integrated autocorrelation function i? Zl i n t(«) for 
P = 1/2, R = 1/4, and po = 3/4 in a ln-ln plot: the upper set 
of points refers to z = p — (p) (solid line is the linear fit), while the 
lower set of points refers to z = (j> — ((f)) (dashed line). 

From our results we draw the following picture of the 
motion in configuration space; cf. Fig. [l]. Starting with 
a small cluster of initial conditions with fixed tangential 
momentum at the boundary away from the supercon- 
ducting interface, the particles move collectively along 
skipping trajectories with constant mean angular veloc- 
ity (vuj(p) divided by the arclength of the orbits) until 
the superconducting region is met. Then, the sign of the 
tangential momentum is inverted. Elementary geome- 
try shows that the mean angular velocity changes. At 
the next collision with the superconducting region the 
same thing happens. Orbits reaching the boundary at 
different sides of a critical point separate after leaving 



the boundary. Hence, the cluster of initial conditions 
starts to spread out on the invariant surface. However, 
weak mixing allows occasional reclusterings with decreas- 
ing frequency (the AF does not decay to zero, but the 
integrated AF does). 

Finally, we address the two limiting cases of zero and 
high magnetic fields. In the zero-field limit, time-reversal 
symmetry is recovered, uj(~- p) — —uj(p). All orbits are 
self-retracing and periodic. The phase space is foliated 
by periodic orbits rather than by two-dimensional sur- 
faces. In the regime of high magnetic field, u>(p) is small 
and close to u(— p). This situation corresponds to the 
symmetric barrier billiard with small p. 



V. CONCLUSION AND OUTLOOK 

We have studied a new kind of pseudointegrable sys- 
tem, the circular Andreev billiard. It is different from 
known pseudointegrable systems in three respects: (i) it 
does not belong to the class of polygonal billiards; (ii) 
it is not symmetric under time reversal; and (iii) it has 
a nontrivial foliation of energy surfaces by invariant sur- 
faces of genus two and two-parameter families of periodic 
orbits (not touching the billiard boundary). We have 
shown that the critical points, i.e., the boundary points 
separating normal and superconducting regions, play the 
role of critical corners in polygons. Moreover, we have 
demonstrated that the dynamics on invariant surfaces 
in the circular Andreev billiard have typically the same 
ergodic properties as in the asymmetric barrier billiard 
(restricted to the chosen Poincare surfaces of section). 
This finding has been used to show that the Poincare 
map is ergodic on the invariant sets. Moreover, we have 
provided numerical evidence that the Poincare map is 
generically weak mixing. 

Weak mixing as maximal ergodic property implies in- 
teresting spectral properties p0| and anomalous trans- 
port. Recent studies in these directions on rational poly- 
gons [ ^p|JTl|Jl^] should be easy to carry over to the circular 
Andreev billiard. 

Of particular interest is the quantum mechanics of 
the circular Andreev billiard because of two reasons. 
First, due to the pseudointegrability we expect an exotic 
quantum-classical correspondence as in rational poly- 
gons |15|. Second, as for rational polygons we might ob- 
serve intermediate energy-level statistics |2l| which are, 
however, modified by the broken time-reversal symmetry 
and the nontrivial foliation. 

The link between the two different classes of systems 
may also throw some light on one aspect of the proximity 
effect in chaotic Andreev billiards, viz., the appearance 
of a gap in the local density of states in an energy in- 
terval above the Fermi energy. Random matrix theory 
can model this gap [^,^3), but the semiclassical the- 
ory predicts an exponential suppression of the density 
of states |24|,^5]]. One origin of this discrepancy could 
be the diagonal approximation used in the semiclassi- 
cal theory E6[. However, our finding indicates another 
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possible explanation. It is well known that the semiclas- 
sical treatment of polygonal billiards requires not only 
periodic orbits but also diffractive orbits, that are orbits 
starting and ending at critical corners; see e.g., p7[ and 
references therein. An interesting research project would 
be to incorporate analogously diffractive orbits (stem- 
ming from the critical points) in the semiclassical theory 
of Andreev billiards and see if this removes the discrep- 
ancy to the random matrix theory. This idea is supported 
by the fact that pointlike scatterers (diffraction) in An- 
dreev billiards gives rise to a gap in the spectrum near 
the Fermi energy pS| ]. 
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